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ABSTRACT 


We report the results of far-infrared and radio molecular CO 
observations of the W3 North star- forming region. The object shows 
extended dust and gas emission, which we interpret as the result of 
fairly advanced disruption of a molecular cloud. This idea is supported 
by an estimate for the age of the embedded H II region. It is possible 
that W3 North is the oldest object among the W3 complex of sources. 

Subject headings : infrared: sources — interstellar: molecules — 

nebulae: H II regions — nebulae: individual (W3) 
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I. INTRODUCTION 


The W3 North (G133.8+1.4) source is the northernmost member of a 

string of active star-forming regions that marks the western boundary of 

the giant H II region W4. Although not often a specific target of 

observational programs, W3N has been included in a number of surveys of 

the region. At far-infrared wavelengths Thronson, Campbell, and 

Hoffmann (1980) found bright, extended emission from the region, 

although they were not able to clearly delineate between internal and 

external heating. Thronson and Price (1982) discussed the large-scale 

mid-infrared structure of the region, within which W3N appeared as a 

shoulder at the edge of the bright W3A source. Lada et al . (1978) 

12 

presented a CO survey of the area at an angular resolution of 216 and 
8‘. These observations clearly place W3N in a ridge of molecular 
emission that includes the three other bright infrared sources, W3A, 
W3(0H), and AFGL. 333. Lada et al . argue that these four sources were 
formed out of material compressed by the expansion of W4. Dickel (1980) 
proposed an alternative model for the source. 

The only infrared observations with an angular resolution 
sufficient to study structure in the source was that of Zeilik (1979), 
although only a 2.2pm map was presented. At this wavelength he found 
fairly weak emission, extended by about l 1 and offset slightly from the 
radio continuum peak. Fazio et al . (1975) scanned close to the location 
of VJ3N and estimated a lower limit to its far-infrared flux density. 
Furniss et al . (1974) mapped much of the W3 complex and detected W3N, 
but used beam sizes of 4 ' -5 ' which made it difficult to analyze the 
object in any detail . 
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Radio continuum surveys of the area (Rohlfs, Braunsfurth, and Hills 
1977; Wendker and Altenhoff 1977) show W3N as a distinct source, but it 
takes angular resolution better than a few arcmin to fully separate the 
many sources in the region (Sullivan and Downes 1973, Schraml and Mezger 
1969). These latter observations suggest that the object has structure 
on a scale size of a few arcmin and is powered by a star — or stars — • 
with a luminosity equivalent to an 07 zero age main sequence (ZAMS) 
star. 

A useful summary of the many observations of the W3/W4 region has 
been presented by Dickel et al . (1980). 

II. THE OBSERVATIONS 
a . ) The Far- infrared Observations 

The observations reported here were obtained at the Nasmyth focus 
of the 0.9m telescope onboard NASA's Kuiper Airborne Observatory during 
January, 1983. The detector system was a multi-filter, seven-element 
array configured with a central detector and six close-packed detectors 
surrounding it. The angular resolution of each of the seven apertures 
was 49" and the center- to-center separation between the middle and an 
outer aperture was 57". The reference beam was 1\ 5 to the northeast of 

W3N. The objects used for calibration were W3(0H) and OMC 1. 

The source was observed in four broad passbands. It was first 

mapped at 60 ym (Fig. 1). The multi-channel array was then centered on 

the maximum emission and we obtained photometry with the three other 
filters. The results are presented in Table 1. The flux density in the 
central beam is presented, along with the flux density measured for all 
seven channels. Only for the 60 ym filter was there sufficient data to 
estimate a total flux density for the source. Inspection of Figure 1 
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suggests that roughly 20% of the total 60 pm flux density was missed by 

our mapping, particularly to the south of the peak. 

The far- infrared luminosity was found for the source by assuming 

that all the filters had the same spatial distribution as that at 60pm. 

No correction was made for possible emission outside the boundaries of 

the area mapped. The object was assumed to be at a distance of 2.3 kpc. 

The flux densities in the central channel at emission peak were 

used to estimate a dust temperature and emission optical depth. 

Following the usual convention for multi-filter far-infrared photometry, 

we assumed that the source spectrum was of the form F « vB (T . ) , where 

v v v d' 

T d is the dust temperature. Although this spectral shape fits the data 
in Table 1, we do not have sufficient spectral resolution and coverage 
to distinguish between the function adopted here or other similar ones. 
[The problem of deriving source parameters from far-infrared photometry 
is considered more fully in Thronson (1985).] Fitting our observations 
by eye gives T^ = 75 ± 15 K, where the uncertainty is our estimate of 
the fit. An emission optical depth at 60pm was found for the peak using 
F v - Q B v (T d )(l - exp[- gg])» where fl is the beam solid angle. For 
the peak of W3N, we calculated = 0.003. 

b . ) The Millimeter-wave Observations 
Millimeter-wave emission from the J = l-*0 transition of CO was 
observed in W3N during September, 1981 using the 11 m telescope of 

r i o 

NRA0. The ^C0 emission was mapped over a 6 ' x6 * area centered on the 

6 The National Radio Astronomy Observatory is operated by Associated 
Universities , Inc. , under contract with the National Science Foundation. 
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13 

far-infrared maximum. A single CO spectrum was obtained at the 

far-infrared maximum. The beamwidth was 66", the map was sampled every 

beam-width, and the reference beam position was 30* to the east in an 

area devoid of CO emission (Lada et al . 1978). The frequency resolution 

for the observations reported here was 250 kHz. The calibration object 

was W3(0H) for which we took T A ( l2 C0) = 16.8 K and T ft ( 13 C0) = 7.1 K. 

The radio molecular data was analyzed in the usual manner (e.g., 

Dickman 1978, Martin and Barrett 1978). The CO emission was assumed 

to be optically thick and its calculated excitation temperature was 

13 13 

assumed also to hold for CO. The CO column density at the 

13 

far-infrared maximum, N( CO), was found under the assumption of 

thermodynamic equilibrium. From inspection of our map of the 
* 12 

distribution of T^( CO), we estimated a forward beam coupling 
coefficient of 0.54. This object is part of an extensive program of 
study of the CO emission from the W3 region. The column density and 
mass distribution are discussed here, but the velocity structure 
analysis is deferred to a later paper (Thronson, Lada, and Hewagama 
1984) . 


III. ANALYSIS AND DISCUSSION 
a . ) Morphology and Energetics 

The luminosity estimated for the source from our far- infrared data, 

4 

7 ± 1.4 x 10 L g , is only slightly below that estimated from radio 
continuum observations (Schraml and Mezger 1969, Sullivan and Downes 
1973). As noted in the previous section, our observations may 
underestimate the luminosity by 20%, so we consider the infrared and 
radio continuum data to be in good agreement. The W3N object is excited 
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by a star — or stars — with a luminosity equivalent to an 07 ZAMS 
star. The agreement between the luminosity determined from the infrared 
and radio observations means that there is very little dust mixed in 
with the ionizing gas and competing for the ionizing photons. This is a 
common situation for evolved diffuse H II regions and is discussed 
further in § I lie. 

The infrared data reported here are also in generally good 
agreement with that found by other investigators, after correction for 
apparent beam-size effects. This means that surrounding the region 
mapped in Figure 1 may be an extended area of low-level emission. In 
addition, with W3A only 12' to the southwest, there can be contamination 
in the large beams of other experiments. In a 4' beam, Furniss et al , 
(1974) estimated a far-infrared luminosity of 1 x 10 L q (all 
observations have been corrected to distance of 2.3 kpc, where 
necessary), in good agreement with our data. However, Thronson, 
Campbell, and Hoffmann (1980) found Lp^ = 1.3x10^ L Q in a 12' beam. 
This latter observation would be particularly susceptible to the bright 
emission from W3A. 

12 

Other than for CO, the maximum emission at all wavelengths 
longward of visual is close to the same position, within the 
uncertainties of observation (2.2 ym and 12.6 urn: Zeilik 1979; 1.4 GHz 
continuum: Sullivan and Downes; 60 ym: this work). There is very little 
information as to size or structure of the source from other 
far-infrared observations. Fazio et. al . (1975) roughly estimated a 
diameter greater than 2', but both Furniss et. al . and Thronson, 
Campbell, and Hoffmann found the object unresolved to their large beams. 
Inspection of Figure 1 shows that 2' is a good estimate of the 
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far-infrared half-power size. This is much larger than the 2.2 urn size 

of the object (Zeilik 1979), an expected result since the near-infrared 

emission is dominated by free-free emission from the ionized gas or by 

the very hottest dust. However, the far-infrared size is only slightly 

larger than that found for the radio continuum (Sullivan and Downes 

* 12 

1973). The distribution of T^( CO) (Figure 2) is much more extended 
than either the far-infrared or radio continuum, which is a fairly 
common situation. It is less common to find only a modest 1 ocal maximum 
at the far-infrared peak, as is the case for W3N. 

Our data, as well as that at other wavelengths, show that VJ3N can 
be described as a partially obscured, diffuse H II region. Visual 
photographs (e.g., Sullivan and Downes; Thronson, Campbell, and 
Hoffmann; Dickel et al . ) clearly give the impression of ionized gas 
breaking through obscuring neutral material. The far-infrared and radio 
continuum observations show that the source has little compact, high 
brightness structure. The extended nature of 60 ym emission suggests 
that the exciting star(s) in the object are not heavily dust enshrouded, 
but instead the stellar and nebular radiation is escaping to heat a 
large area of surrounding material. The fairly low molecular gas 
temperatures (Table 1), even at the far-infrared peak, is further 
evidence for more extended, diffuse structure in the source. 

b ) Mass and Density 

Because the dust temperature is so high (Table 1) and the flux 
density is only moderately high, the 60 ym emission optical is low. 
This supports our description of W3N as a diffuse H II region in which 
the typical far-infrared-emitting grain is embedded in a more intense 
radiation field than would be the case for a dense, compact H II region. 
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The low emission optical depth is reflected in the visual appearance of 

13 13 

the source, as well as the modest CO column density, N( CO) = 7 x 

15 

10 ~ cm . Using the relations in Dickman (1978), this is equivalent to 
3 visual magnitudes of extinction. This level of obscuration, although 
low, is sufficient so that most of the visual and ultraviolet photons 
from the source will be absorbed by dust and re-radiated at infrared 
wavelengths. This is necessary to explain the agreement between the 
luminosity determined from the radio continuum observations and from the 
infrared data. 

As Figure 2 shows, the molecular cloud associated directly with W3N 

is not clearly separated from more extensive material in the area. 

Therefore, only a rough estimate of the mass of the molecular material 

specifically associated with the infrared source can be made. If we 

take the W3N cloud to be 3' (2 pc) in diameter, the relationship between 
13 

mass and CO column density found by Dickman results in a cloud mass of 

13 

230 M . This assumes that the CO column density at the infrared peak 

is characteristic of the rest of the source as well, which is unlikely. 

It is generally assumed that in most molecular clouds near active 

star-forming regions, the gas is heated by collisions with dust. We 

therefore expect that the material that dominates the far- infrared 

emission is mixed with that dominating the CO emission. 

13 

The 60 y m optical depth should be related to the CO column density. 

-19 13 

From our data at the infrared maximum we find = 4 x 10 N( CO) = 
0.001 A . This result is in good agreement with similar calculations by 
a number of investigators (see discussion in Thronson et al . 1984). 

A fundamental parameter for any star- forming region is the volume 
density and we use the column density to calculate it for W3N. If the 
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1o r 

object is 2 pc in diameter and we use N(H 2 )/N( CO) = 5 x 10 from 

a- _3 

Dickman, we calculated n^ ^ 600 cm from the CO column density. This 
fairly low level is consistent with our picture of the source as a 
diffuse H II region. 

c . ) The Age of W3 North 

It is the visual appearance of the object (Lada jet _aj_; Thronson, 

Campbell, and Hoffmann) that most quickly gives the impression that this 

young, star-forming region has substantially disrupted its surrounding 

neutral cloud. The object shows extensive, visually bright emission 

centered very close to the infrared and radio continuum maximum. The 

12 

fact that the 60 y m and CO emission are both very extended further 
supports the idea that the cloud is being dissipated. 

As part of an extensive study of the evolution of the W3 system 
(see Thronson, Lada, and Hewagama for a more complete discussion), it is 
valuable to estimate an age for W3N. There is really only one technique 
suitable for this: deriving the expansion age of the H II region. This 
requires knowing initial conditions in W3N , but the uncertain values 
enter to a low power in the final expression. Spitzer (1978) derived 
equations describing the expansion of an H II region, from which we may 
write its age 


L a, 


4 r 

Tv 


7/4 



Here r is the present radius of the ionized zone (0.6 pc, Harten 1976), 
V is the expansion velocity (^10 km s~^), and r g is the Stromgren radius 
of the H II when expansion began (with r g «r). This latter parameter 
may be replaced, using the standard expression for the Stromgren sphere. 
We therefore get - - - 


t % 7 


4Tra 
3TT , 

i 


,7/4 n 1/2 
V 
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where a is the hydrogen recombination coefficient, is the Lyman 

continuum flux for an 07 star, and n is the electron density when the H 

e 

II region was originally established. It is clear that n g is the most 

poorly known parameter, but on the basis of the most compact H II 

5 -3 

regions discussed in the literature, we adopt n % 10 cm . The 

5 

resulting estimate for the age is t > 4 x 10 years. We consider this a 
weak lower limit since, of the parameters in the above equation most 
subject to change, V may be lower than adopted here and n g may be 
1 arger . 

Despite the uncertainties in t, it is worth noting that our 
estimate of the age of the H II region is very close to the age that 
Dickel (1980) estimates for the time that higher-density "clumps" (not 
stars) existed in the W3 region. In other words, the estimate for the H 
II region that we derive here suggests that at least some stars in the 
complex were more evolved than in the Dickel scenario. 
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IV. SUMMARY 


The luminous W3 North H II region has been mapped at far- infrared 

12 

wavelengths and in the emission of the J = l->0 line of CO. We deduce 
the following characteristics about the source. 

1. The object is excited by the equivalent of an 07 ZAMS star and 
is surrounded by extended CO and far-infrared emission. 

2. Only a moderate-mass neutral cloud surrounds the H II region. 

_3 

Its density is in the range of n H ^ 600 cm and appears to have an 

2 

average visual extinction of <u3 mags. 

5 

3. The expansion age of the H II region is t ^ 4 x 10 years. 
This is long enough to significantly disrupt the overlying material, 
which is exactly how the source appears in visual photographs. 
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TABLE 1 


SOURCE PARAMETERS 


Peak Position 


a(1950)=2 h 23 m 

1 ?5±2 S 

(60 m) 


6 (1950)=62° 2' 

10 "±15" 

Flux Density a (Jy) 

Peak 

7-Channel s 

Total 


(49" Beam) 


40ym 

1100 

3600 

- 

60ym 

910 

3600 

4900 

100 urn 

730 

3200 

- 

160ym 

250 

1200 

- 

Luminosity (30ym - 200ym)^ 


7 x 10 4 L q 


Peak Dust Temperature 


T , = 75±15 K 
d 


Peak Gas Temperature 


T* 1 ? 

'A( ^CO) = 12.5 

K 



T A( l3 C0) = 1.5 

T ex " 29 K 

K 

Peak 60 m Optical Depth 


0.003 


13 

Peak CO Column Density 


7 x 10 15 cm" 2 


Central Volume Density 0 


600 cm" 3 


Cloud Mass 


230 M 

0 


a 

Internal uncertainty less than 

1%; Calibration uncertainty ± 

20%. 

b May be an underestimate by about 

20%; see § 

I la . 



Determined as described in § lib. 
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FIGURE CAPTIONS 


Fig. 1 - The distribution of 60 ym flux density in W3 North. The 

contours are normalized to 1.00 = 910 Jy. The light outer boundary 
shows the extent of the area mapped. The beam size is represented 
by the cross-hatched circle. The cross is located at the emission 
maximum. The length of the four arms are each 15" long, the 2a 
uncertainty in absolute position. 

* 12 

Fig. 2 - The distribution of J = 1 -*■ 0 T^ CO. The contours are labeled 
in kelvins. The beam size is represented by the cross-hatched circle. 
The cross is located at the position of maximum 60 ym emission. The 
arms of the cross are each 15" in length, the pointing uncertainty 
of the 11 m telescope. 
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